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Abstract 

The instantaneous contributions to the third post-Newtonian (3PN) gravitational wave luminosity from 
the inspiral phase of a binary system of compact objects moving in a quasi elliptical orbit is computed us- 
ing the multipolar post-Minkowskian wave generation formalism. The necessary inputs for this calculation 
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orbits the flux is averaged over the binary's orbit. Supplementing this by the important hereditary contri- 
butions arising from tails, tails-of-tails and tails squared terms calculated in a previous paper, the complete 
3PN energy flux is obtained. The final result presented in this paper would be needed for the construction 
of ready-to-use templates for binaries moving on non-circular orbits, a plausible class of sources not only 
for the space based detectors like LISA but also for the ground based ones. 
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I. INTRODUCTION 



Inspiralling compact binaries, one of the prototype sources for laser interferometric gravita- 
tional wave (GW) detectors, are usually modelled as moving in quasi-circular orbits. This is 
justified since gravitational radiation reaction, under which it inspirals, circularizes the orbit to- 
wards the late stages of inspiral |Q1 @]. This late phase of inspiral and the ensuing merger phase 
offer promises for the GW interferometric detectors. The recently discovered double pulsar sys- 
tem (j, 3] has an eccentricity as low as 0.088 consistent with the circular orbit assumption for 
the late inspiral and pre-merger phases, believed to be reasonable enough for most of the binary 
systems made of neutron stars or black holes (BHs). 

The theoretical modelling of the binary's phase evolution to a very high precision is called 
the phasing formula. This is the basic theoretical ingredient used in the construction of search 
templates for GW using matched filtering (H]. The two key inputs required for the construction 
of templates for binaries moving in quasi-circular orbits in the adiabatic approximation are the 
orbital energy and the GW luminosity (energy flux). These are computed using a cocktail of 
approximation schemes in general relativity. The schemes include the multipole decomposition, 
the post-Minkowskian expansion of the gravitational field or non-linearity expansion in Newton's 
constant G, the post-Newtonian expansion in v/c, and the far-zone expansion in powers of \/R, 
where R is the distance from the source (see [6] for a recent review). 

Though the garden variety binary sources of GWs for terrestrial laser interferometric GW de- 
tectors are those moving in quasi-circular orbits, there is an increased recent interest in inspiralling 
binaries moving in quasi-eccentric orbits. Astrophysical scenarios currently exist which lead to bi- 
naries with non-zero eccentricity in the GW detector bandwidth, both terrestrial and space-based. 
For instance, inner binaries of hierarchical triplets undergoing Kozai oscillations f7] could not 
only merge due to gravitational radiation reaction but a good fraction (~ 30%) of them will have 
eccentricity greater than about 0.1 as they enter the sensitivity band of advanced ground based 
interferometers Almost all the above systems possess eccentricities below 0.2 at 40 Hz and 
below 0.02 at 200 Hz. The population of stellar mass binaries in globular clusters is expected to 
have a thermal distribution of eccentricities [@]. In a study on the growth of intermediate BHs ifioll 
in globular clusters it was found that the binaries have eccentricities between 0.1 and 0.2 in the 
LISA bandwidth. Though, supermassive black hole binaries are powerful GW sources for LISA, 
it is not yet conclusive if they would be in quasi-circular or quasi-eccentric orbits 111 ill . If a Kozai 
mechanism is at work, these supermassive BH binaries could be in highly eccentric orbits and 
merge within the Hubble time II 1211 . Sources of the kind discussed above provide the prime moti- 
vation to investigate higher post-Newtonian order modelling for quasi-eccentric binaries. 

The GW energy flux or luminosity from a system of twopoint masses in elliptic motion was 
first computed by Peters and Mathews at Newtonian order njj, |2Q . The post-Newtonian (PN) cor- 
rections to the gravitational wave flux at 1PN and 1.5PN were provided in El El EE [3 and 
used to study the associated evolution of orbital elements u sing the 1PN "quasi-Keplerian" repre- 
sentation of the binary's orbit lfl8ll . Gopakumar and Iyer [U 20j] further extended these results 
to 2PN order using the generalized quasi-Keplerian representation developed in Ref. 112 ll l22l 12311. 
The results for the energy flux and waveform presented in |l9l]_was in perfect agreement with those 
obtained by Will and Wiseman using a different formalism II24I1 . Recently, Damour, Gopakumar 
and Iyer [25] discussed an analytic method for constructing high accuracy templates for the GW 
signals from the inspiral phase of compact binaries moving on quasi-elliptical orbits. They used an 
improved "method of variation of constants" to combine the three time scales involved in the el- 
liptical orbit case, namely, orbital period, periastron precession and radiation reaction time scales, 



2 



without making the usual approximation of treating the radiative time scale as an adiabatic process. 

The generation problem for gravitational waves at any PN order requires the solution to two 
independent problems. The first relates to the equation of motion of the binary and the second 
to the far zone fluxes of energy, angular momentum and linear momentum. The latter requires 
the computation of the relativistic mass and current multipole moments to appropriate PN orders. 
The 3PN equations of motion (EOM) required to handle gravitational wave phasing turned out to 
be technically very involved due to the issues related to the self-field regularization using Riesz 
or Hadamard regularizations [26, [27b - Only by a deeper understanding of the origin of the am- 
biguities in Hadamard regularization, and the use of dimensional regularization has the problem 
been uniquely resolved 11281 . 12911 and provided the value of the ambiguity parameter oj s n26\ or 
equivalently A I127ll . We thus have in hand the requisite 3PN EOM for compact binaries moving in 
general orbits. The computation of the GW luminosity at 3PN or (v/c) 6 beyond the leading Ein- 
stein quadrupole formula crucially requires the computation of the 3PN accurate mass quadrupole 
moment. For its completion the same technique as in the EOM was successfully applied, namely, 
to compute using Hadamard's regularization all the terms except a few terms parametrized by am- 
biguity parameters (which turn out to be three, denoted k and £) [30, 3 ID, and then to determine 
the value of these parameters by computing the difference between the dimensional and Hadamard 
regularizations [31, |3_2j,|3_3j,|3_4D. These works thus provide the fully determined 3PN accurate mass 
quadrupole for general orbits - the other important ingredient to compute the 3PN accurate energy 
and angular momentum fluxes for inspiralling compact binaries moving in general non-circular 
orbits. The 3.5PN phasing of inspiralling compact binaries moving in quasi-circular orbits is now 
complete and available for use in GW data analysis 11321 13511 . Note that the 3PN contribution to 
the energy flux comes not only from the "instantaneous" terms discussed in this paper but also 
includes "hereditary" contributions arising from tails, tails of tails and tail- squared terms. A semi- 
analytical scheme is proposed and discussed in detail in a companion paper II36I1 1 to evaluate these 
history-dependent contributions. 

In this paper, for binaries moving in elliptical orbits, we compute all the instantaneous contribu- 
tions to the 3PN accurate GW energy flux. The orbital average of this flux will be obtained using 
the 3PN quasi-Keplerian parametrization of the binary's orbital motion recently constructed by 
Memmesheimer, Gopakumar and Schafer [371]. We shall supplement these by contributions from 
the hereditary terms computed in Paper I. The final expression will represent gravitational waves 
from a binary evolving negligibly under gravitational radiation reaction, including precisely up to 
3PN order the effects of eccentricity and periastron precession during epochs of inspiral when the 
orbital parameters are essentially constant over a few orbital revolutions. It also represents the first 
step towards the discussion of the qua si- elliptical case: the evolution of the binary in an elliptical 
orbit under gravitational radiation reaction. The present work extends the circular orbit results 
at 2.5PN 113 811 and 3PN DO. 13211 to the elliptical orbit case. Further, it extends earlier works on 
instantaneous contributions for binaries moving in elliptical orbits at 1PN Ill4lll5ll and 2PN IU9I1 
to 3PN order. Similarly, Paper I extends hereditary contributions at 1.5PN lllql to 2.5PN order 
and 3PN, where the 3PN hereditary contributions comprise the tails of tails and are extensions of 
Refs. ff39L l40Tl for circular orbits to the elliptical orbit case. 

In Sect. [II] we begin with the structure of the far-zone flux of energy, use expressions relating 
the radiative moments to the source moments and decompose the energy flux expression into 
its instantaneous and hereditary parts. Section [HI] lists all the requisite multipole moments in 
standard harmonic coordinates for binaries moving in general (non-circular) orbits. Section [IV] 



1 Hereafter Ref. 136(1 will be called Paper I. 
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introduces the 3PN equations of motion which are necessary to handle the time derivatives of the 
moments. Section |V] discusses the computation of the instantaneous terms in the energy flux and 
Sect. [VI] recasts the flux in modified harmonic coordinates (without logarithms at 3PN order) and 
Arnowitt, Deser and Misner (ADM) coordinates. Section |VTH summarises the 3PN quasi-Keplerian 
representation required to average the flux expression over an orbit. Section IVIIII exhibits the 
orbital average of the energy flux in modified harmonic coordinates and ADM coordinates, and 
finally provides an expression of the complete energy flux in terms of gauge invariant variables. 



H. THE FAR ZONE FLUX OF ENERGY 



In this section, we discuss the computation of the 3PN accurate energy flux for general isolated 
sources. Starting from the expression for the far zone flux in terms of the radiative multipole mo- 
ments and using the relations connecting the radiative multipole moments to the source moments, 
we write the resultant structure of the GW energy flux. 

Following Thorne 114 ill , the expression for the 3PN accurate far zone energy flux T = 
(d£/dt) GW in terms of symmetric trace-free (STF) radiative multipole moments reads as 2 
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In the above Ul and Vl (where L = i\i 2 ■ ■ ■ U represents a multi-index composed of / spatial indices) 
are the mass-type and current-type radiative multipole moments respectively and U® and V® 
denote their I th time derivatives. The moments are functions of retarded time U = T - R/c in 
radiative coordinates. 

In the multipolar-post-Minkowskian (MPM) formalism, the radiative moments U L and V L can 
be re-expressed in terms of the source moments to an accuracy sufficient for the computation of 
the energy flux. For the flux to be complete up to 3PN approximation, one must compute the mass 
type radiative quadrupole £/ ;y to 3PN accuracy, mass octupole and current quadrupole Vjj to 
2PN accuracy, mass hexadecapole f//^ m and current octupole to 1PN accuracy and finally 
Uijkmn an d Vijkm to Newtonian accuracy. 

The relations connecting the different radiative moments Ul and Vl to the correspon ding so urce 
moments I L and J L are given below. For the 3PN mass quadrupole moment we have B8U39lI401 |42|1 
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where the bracket <> surrounding indices denotes the STF projection, and e a u is the usual Levi- 
Civita symbol such that e 12 3 = +1 . The 7 L 's and J L 's are the mass and current-type source moments 
(and if , denote their p th time derivatives), and W is the monopole corresponding to the set of 
"gauge" moments Wl, using the same definitions as in lf30h . In the above formula, M (which is in 
factor of the tail integral at 1.5PN order and the tail-of-tail integral at 3PN) is the total ADM mass 
of the source. The non-linear memory integral at 2.5PN is a time anti-derivative and will become 
instantaneous in the energy flux. The moments needed at 2PN order include only the dominant 
tails and are 
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For all the other moments required in the computation we need only the leading order accuracy in 
the relation between radiative and source moments, so that 



U L (U) = I®(U)+0(3), 
V L (U) = J ( ?(U) + 0(3) . 
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The constant length r scaling the logarithm is the one introduced in the general MPM for- 
malism and has been chosen here to match with the choice made in the computation of tails-of- 
tails in ll40ll . It is a freely specifiable constant, entering the relation between the retarded time 
U = T - R/c in radiative coordinates and the corresponding retarded time tu - r^/c in harmonic 
coordinates (where r H is the distance of the source in harmonic coordinates). More precisely we 
have 



c 



2CM l„^i+0(5). 



(2.5) 



From Eqs <12.2D - (12.4D , it is clear that the radiative moments have two distinct contributions. One 
part depends on the moments only at the retarded time, U = T - R/c; this part is referred to as the 
"instantaneous contribution", and forms the subject matter of the present paper. The second part 
on the other hand depends on the dynamics of the system in its entire past, i.e. at any U - r < U, 
and is referred to as the "hereditary contribution". Equally important but requiring a different 
treatment, the hereditary contribution is dealt with in Paper I as mentioned earlier. We are thus 
allowed to write down explicitly the different kinds of contributions to the far zone energy flux up 
to 3PN. We have, 

J = Tmst + ^hered > (2.6) 

where the instantaneous contribution of interest in this paper is explicitly given by 
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The hereditary contribution is given in Sect. Ill A of Paper I. We recall that it is decomposed as 
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All the tail contributions are thoroughly computed in Paper I and we shall use those results to 
obtain the complete GW energy flux in Sec. IVIIII 



III. THE MULTIPOLE MOMENTS OF COMPACT BINARY SYSTEMS 



We provide, in this Section, the requisite multipole moments needed for the computation of 
the 3PN accurate energy flux for compact binaries in the standard harmonic coordinate system. 
By standard harmonic coordinates we refer to the specific coordinate system which has been used 
consistently in previous works 11271 |29j, m sums |35|, |43|]. We recall that these coordinates 
contain some logarithms at the 3PN level both in the equations of motion of the binary 11271 . 14311 
and in their multipole moments l!30Ll3lL 13211. Later, we shall also define some modified harmonic 
coordinates which do not involve such logarithms at the 3PN order. 

The multipole moments are generalisations to non-circular orbits of the expressions available in 
Ref. [30] for circular orbits. They are computed by implementing the detailed method described in 
Re f. 113 ill. Though algebraically long and involved, the procedure is fairly algorithmic as explained 
in [ 3(1" 31]. We thus skip all those details of computations and list the final results we need. The 
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3PN mass quadrupole /y is already given in Ref. [3V\\ and its expression (valid in the frame of the 
center of mass) is 
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In the above equation r is the length scale appearing in the definition of the source multipole mo- 
ments 113 ill and is the same as in Eq. ( 12.51) . On the other hand, the different constant r' Q is related to 
two other length scales r\ and r' 2 (one for each particle) by m In r' = m x In r[ + mi In r' v and is spe- 
cific to the application of the formalism to point particle systems. It comes from regularizing the 
self-field of point particles in the standard harmonic coordinate system. It is very important to note 
that the two length scales r[ and r' 2 are the same as the two scales that appear in the final expression 
of the 3PN equations of motion in standard harmonic coordinates 12711 . The requirement that these 
r'j and r' 2 should match with similar scales that appear in the equations of motion determine, using 
dimensional regularization, the values of the Hadamard's regularization constants k and f that 
formerly appeared in the 3PN multipole moments |3(1 31]. The regularization constants are thus 
determined and we have consistently replaced k and f by their values known from 11321 . 13411 . The 



constants r[, r' 2 and hence r' Q are "unphysical" in the sense that they can be arbitrarily changed by 
a coordinate transformation of the "bulk" metric outside the particles Q27Q, or, more appropriately 
(when considering the renormalisation which follows the dimensional regularization), by some 
shifts of the particles' world lines |29[ 34ll . 

The 2PN mass octupole and current quadrupole moments for general orbits are the other non- 
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trivial moments required. They are given by 
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In the above and what follows, Xjj k ... = XiXjX k • ■ • and Vjj k ... = ViVjV k 
the STF projection. The 1PN moments read as 
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[The underlined index a means that it should be excluded from the STF projection.] Finally we 
also need 
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as well as W, the monopole corresponding to the gauge moments W L , and which is given by 
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W = -vmrr + 0(2). 
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(3.6) 



IV. THE EQUATIONS OF MOTION OF COMPACT BINARY SYSTEMS 
A. The equations of motion in standard harmonic coordinates 



The computation of the flux will involve the time derivatives of the latter source moments. 
The 3PN accurate flux requires the 3PN equations of motion for compact binaries which are now 
complete 11271 . |28L |29l l44l 14511 . For the present work, where the multipole moments are computed 
in standard harmonic coordinates and reduced to the centre of mass (CM) frame, we require the 
3PN accurate equation of motion (or acceleration) in the CM frame associated with the standard 
harmonic gauge. This was computed in Ii43ll and given as 



a i = ^ = -^[(l + P)n i + Qv i \ + 0(7), 



where the coefficients P and Q are: 



(4.1) 



P=^<- 



1 f 3r 2 v 



Gm 



+ i/ + 3v(/ (4 + 2v) 

2 r 



1 (I5r 4 v 45r 4 v 2 9r 2 vv 2 



+ 7 



+ 



+ 



8 

Gm 
r 

G 2 m 2 



8 2 
2 r 2 - 25 r 2 v - 2 r 2 v 
87 v 



+ 6 r 2 v 2 v 2 + 3 v v 4 - 4 v 2 v 4 
13 vv 2 



+ 2v 2 v 2 



4 



1 f 24rvv 2 Gm 136rvG 2 m 2 

+ c E \ 5^ r 15 ~ 

If 35r 6 v 175 r 6 v 2 175 r 6 v 3 15 rW 
+ — + — — — + 



16 



16 



16 



135 r 4 v 2 v 2 255r 4 vV 15 rW 237 r 2 v 2 v 4 



+ 



4 8 2 

45 r 2 vV 11 vv 6 49v 2 v 6 



2 
+ 



+ 



Gm / , 
+ — 79 r 4 v - 
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69r 4 v 2 



4 



+ 13 v 3 v 6 



-30rV- 121 r 2 vv 2 + 16r z v 2 v 



2,., .2 



JJ „2 
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+ 2 0r 2 vV + -^^ + 8vV-10vV 



G 2 m 2 



+ 



4 

, 32573 r 2 v llr 2 v 2 , - 615 r 2 vn 2 26987 vv 2 
r + + 7 r v + ■ 



+ v 3 f; 2 



168 8 

123 V7T 2 v 2 



64 



110r 2 v ln(-) + 22vv 2 ln(-) 



64 

r 

7 

' n 



G 3 m 3 



+ 



437 v 71v 2 41 vtt 2 
■16 - — + 



2 = 4{-4r + 2ry} 



c 

1 (9r 3 v 



+ — 



4 



15 rvo 2 



16 



+ 



2 
Gm 



41/ " V a 7 

2r + + 4 r v 



+ 7 



1 (Svv 2 Gm 24 vGV 
— + — 



5 r 
45 r 5 v 



8 



+ 15 r 5 v 2 + 



15r 5 v 3 



+ 12rW 



111 f 3 v 2 v 2 ,^^^9 65 fvv 4 „ „ 94 , , 4 
■ 12r 3 v 3 v 2 + 19rv 2 y 4 + 6rv 3 i; 4 



+ ■ 



Gml 329 r 3 v 59r 3 v 2 



+ 



8 



+ 18rV - 15rvy 2 -27rvV- 10rv 3 v 2 



G 2 m 2 



+ 



18169rv _ 2 n , 123rv^- 2 
-4 r — — + 25 r v l + 8 r v - 



840 



32 



+44 rv ln(— ) 



(4.2a) 



(4.2b) 



Recall that there was initially a regularization ambiguity constant denoted A in |27|], which has 



been replaced here by its uniquely determined value A = -^§5 |29j]. On the other hand the 
constant r' Q is the same as the one in the 3PN quadrupole moment (I3.1I) - (I3.2I) . 



B. The modified harmonic coordinates (without logarithms) 

The standard harmonic (hereafter SH) coordinate system used till now is useful for analytical 
algebraic checks, but contains gauge-dependent logarithmic terms that are not very convenient in 
numerical calculations. More importantly, in the presence of the logarithmic terms the simple 
generalized quasi-Keplerian representation (reviewed in Sect. IVIII) is not possible impeding the 
process of averaging the flux over the orbital period. Consequently, it is useful to have the expres- 
sion for the energy flux in a modified harmonic (MH) coordinate system without logarithms like 



the one explicitly used in 114611 (we shall alternatively use ADM type coordinates which are also 
free of such logarithms at 3PN). This will require us to re-express the instantaneous expressions 
for the energy flux [given by Eqs. (15.21) below] in terms of corresponding variables in the MH or 
ADM coordinate systems. We provide in this section the definition of the MH coordinate system. 
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Consider the coordinate transformation x' M = x^ + e^(x) which removes the logarithms \n(r/r' ) 
at the level of the equations of motion as discussed in Ref. Il27ll . It is given by 



22 G 2 m,[ m 2 



(4.3) 



where r\ = |x - yj and r 2 = |x - y 2 | are the distances to the two particles with trajectories y\(t) 
and y' 2 (t), and where r = |yi - y 2 | is their relative distance. Following ll29ll the logarithms can be 
equivalently removed by the shifts (sometimes also called the "contact" transformations) of the 
particle world-lines induced by the change of coordinates, namely 



(4.4a) 
(4.4b) 



The equalities here are functional relations, i.e. the two sides of the equations are evaluated at 
the same coordinate time say, t. The spatial shifts £j and £ 2 of the two world-lines are related to 
the coordinate transformation restricted to the world-lines [denoted d[(t) = Ef(t, y\(t)) and e^(t) = 
^(t,y 2 i.t))} by 



^ = £ ^ v l £ o +0 ( s 2 ), 



(4.5a) 
(4.5b) 



where v\ = dy\/dt and v' 2 = dy' 2 /dt are the coordinate velocities. The latter relations are valid 
at linear order in ef[ and di. Now the coordinate transformation (14.31 ) is at 3PN order so we have 
e° = 0(7) and s' = 0(6). Hence, we see from (14.51 ) that at the 3PN order the shifts simply agree 
with the spatial components of the coordinate transformation, 



£=4 + 0(8) 



■0(8) . 



They are readily obtained from Eq. (14.31) as 



22 G 3 m 2 m 2 



c 6 r 2 \ r 



n> In - +0(8), 



22 G 3 mi m 2 



rrin, I r . 



(4.6a) 
(4.6b) 



(4.7a) 
(4.7b) 



Under the shifts of world-lines the accelerations of the particles are changed by the amounts 
8^a\ and 6^a' 2 (i.e. such that the functional equalities a'\ 2 = a\ 2 + ^a'j 2 hold) given by 



8^a\ 
s C a 2 



dt 2 
d% 

dt 2 



-(4-«H(^)+<w>, 



+ 



G/ni 



+ 0(8) 



(4.8a) 
(4.8b) 
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where the second terms come from re-expressing the gravitational force - gradient of the New- 



tonian potential - in terms of the new trajectories (14.41) . The relative acceleration a 1 
changed by the amount 



a\ - a\ is 



d%2 
dt 2 



Gm 



0(8), 



(4.9) 



where m = m x + m 2 and g\ 2 = g\ 
acceleration associated with the shifts ( 14.71) is 



£L 3 An easy calculation shows that the change in the relative 



6^a' = 



G 3 m 3 
c 6 r4 



(-11 Or 2 + 22v 2 )n i + 44rv> 



+ 



176 , 22 , 22 Gm 

r v + 

3 3 3 r 



44 

n" rv 

3 



— riA +0(8) 



(4.10) 



Adding the above shift to the expression for the relative acceleration in SH coordinates as given by 
Eqs. (I4.1I) - (I4.2I) . yields the expression for the acceleration in MH coordinates. Since a" = a 1 + S^a' 
is a functional identity, the resulting MH acceleration is obtained as a function of the "dummy" 
variables denoted v 2 , r and r. Evidently these variables are to be interpreted as the natural variables 
describing the binary motion in MH coordinates. 4 As expected, the logarithms in Eq. (14.101 ) exactly 
cancel the logarithms in the SH acceleration (I4.1I) - (I4.2I) . Some 3PN coefficients in the EOM are 
also modified and the final result agrees with that displayed in Ref. [46] (see also 



For completeness we note also that the above shifts will modify the 3PN conserved energy of 
the binary (associated with the conservative part of the 3PN equations of motion) by the amount 



For the case at hand with the shifts (14.71) and in the center-of-mass frame we find 



(4.11) 



22 G 3 m 4 v 2 
1 c 6 r 3 



Gm 



-3r 2 + v 2 



In - +r 2 +0(8) 



(4.12) 



Comparing with the 3PN energy in SH coordinates as given by Eq. (4.8) in |43|], we see that 
the logarithms ln(r/rg) are also cancelled in the expression for the energy by going to the MH 
coordinates. 



V. THE INSTANTANEOUS PART OF THE 3PN ENERGY FLUX 

Using the multipole moments given in Eqs. (|3.1| )- (I3.5I) , one computes the required time deriva- 
tives with the help of the equations of motion (14.11 ) and obtains the instantaneous part of the energy 
flux as defined by Eq. (I2.7I ). Here we are working in SH coordinates, in which the equations of 
motion are given by Eqs. (I4.1I) - (I4.2I) . In the next section we consider the case of alternative coordi- 
nate systems. The hereditary part computed in Paper I will be added after the process of averaging 



3 This means that x' MH = x' SH + %' n . 

4 To avoid making the notation too heavy we do not add a subscript MH on the variables u 2 , r and r. In the following 
our notation may not always be completely consistent but should be clear from the context. 
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over one orbit (this contribution is the same in all alternative coordinate systems considered in this 
paper). Though lengthy the computation of the different parts constituting the instantaneous terms 
in the energy flux at 3PN order is straightforward. 5 We write the result as 



r** = rz t + rg* + r£« + r^r + + oo) . (5.D 

and find that the various PN pieces are given by 



mst 5 c 5 r 4 { 12 
, PN 32 G 3 m 4 v 2 ( 4 /785 71 \ ,,/ 1487 58 
-r = y -^r Y (336 " 28 V ) + ' V ("W + T 

Gm 2 ( 170 10 \ 4 /687 155 
+ it — — + —v + r -— - — — v 



r \ 21 21 / V 1 12 28 
Gm 2 /367 5 \ G 2 m 2 ( 1 4 

2PN _ 32 G'»V f J47 5497 2215 ., 



inst 5 c 9 r 4 i \14 504 252 

573 1713 1573 , 
+rc — — + — — - v - 



56 28 42 
Gm A 247 5237 199 , 

+ V + V V 



+/•' 



14 252 36 

4 2 / 1009 5069 631 , 

rv v + v 

\ 84 56 14 

Gm 2 2 /4987 8513 2165 , 

+ r v —— — — v + — — v 

r \ 84 84 84 

G 2 m 2 2 /281473 2273 13 , 

+ — V 1^072" + ^5Y V+ 27 y 

2501 10117 2101 7 

+r° + v v 

504 252 126 

Gm A 5585 60971 7145 , 
+ r — + _., v - ——v 



r \ 126 756 378 

G 2 m 2 J 106319 1633 16 , 

H r v v 

r 2 \ 3024 504 9 

G 3 m 3 / 253 19 4 2 \\ 

+ — —\ + —v y 2 }}, (5.2c) 

r 3 \ 378 7 27 // v 

_2spn 32 G 3 m 4 v 2 f / 12349 Gm 4 4524 Gm , , 2753G 2 m 2 2 

^rz.sra _ i j. i 4 + / 

5 c 10 r 4 I I 210 r 35 r 126 r 2 



5 In order to perform some independent checks on the long and involved algebra, we have found it expeditious to 
make two computations using the two harmonic coordinate systems: SH containing the (gauge dependent) log 
terms a la I43I1 and MH without log terms as in Refs. 
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The new results are the instantaneous terms at 2.5PN and 3PN orders. Up to 2PN order, all the 
terms match with those obtained in Refs. [Q], \\4, 19]. As one may notice, the 2.5PN terms in the 



above equation are all proportional to f and hence are zero for the circular orbit case in agreement 
with the result of |[38|l . The r dependence of the 2.5PN terms is important when we discuss their 
orbital average in Sect. I Villi The 3PN terms provide the generalization of the circular orbit results 
in Ref. ll3~oll . As expected, the constant r present in the expression of the mass quadrupole moment 
appears in the final expression for the 3PN flux (the presence of r' Q is of a different type and is dealt 
with in the next Section). The dependence of the instantaneous terms on the scale r should exactly 
cancel a similar contribution coming from the tail terms as determined in Paper I. This cancellation 



has already been checked for circular orbits in 1130m and we shall prove this cancellation for quasi- 
elliptical orbits in Sect. IVIII1 



VI. THE 3PN ENERGY FLUX IN ALTERNATIVE COORDINATES 

The dependence on r' of the result (15.1I) - (15.2I) is due to our use of the SH coordinate system. For 
circular orbits, it was shown [30] that this r' Q dependence disappears when the total flux is expressed 
in terms of the gauge invariant parameter x = (Gmto/c 3 ) 2/3 related to the GW frequency. In the 
general orbit case we shall transform away the dependence on r' by going to different coordinate 
systems such as the MH coordinates studied in Sect. IIVB1 Subsequently we shall average the 
energy flux over an orbital period and exhibit alternative representations of the energy flux for 
elliptical orbits. In particular some of these are in terms of gauge invariant variables related to 



those suggested in Ref. 113711 . 



A. The modified harmonic coordinates 



We now provide the energy flux T in the modified harmonic (MH) system avoiding the appear- 
ance of the logarithms \n(r/r' ) at 3PN order and which has been introduced in Sect. IIVB1 First 
we notice that T is a function of the "natural" variables r, r and v 2 [see Eqs. (15.21) 1. and is a scalar, 



therefore it satisfies, under the shifts of these variables defined by (14.41) . 



T[r,r,v 2 ] = T' W,r', v' 2 ] . 



(6.1) 



This means that we shall have the functional equality T' =< F + 8(F in which 



(6.2) 



where 



S ( r = n^\ 
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Jl2 v -rn 
6 f r = ri + £ 



6 f v 2 



2v 



dt 

i J 12 

dt 



12 



\+o{e). 



(6.3) 



16 



(Recall £' 12 = g\ - g' 2 .) Since the previous formulas are at linear order in the shifts and we are 
interested in the 3PN approximation, they are valid for any shifts at the 3PN order (case of the 
MH coordinates) and also at the 2PN order like the ones associated with the passage to ADM 
coordinates - in the latter case, the error will be at 4PN order. 

In the case of the MH shifts, which start at the 3PN order, one can make an alternative computa- 
tion of the modification of the energy flux. Indeed the only modification vis a vis the calculation in 
standard harmonic (SH) coordinates is the one related to the mass quadrupole moment which must 
be computed to 3PN accuracy. Under the shifts of the particles' trajectories y\ 2 as given by (14.41) . 

the mass quadrupole moment which equals Ijj = mi y^' J> + m 2 y 2 ' ]> + 0(2) in the Newtonian 
approximation, is shifted by the amount (1 «-» 2 meaning the same term but for the other particle) 



2mi yf^f + 1^2 + (9(8) . 



(6.4) 



where the remainder (9(8) comes from the 1PN corrections in the quadrupole moment coupled with 
the 3PN shifts. Using the explicit expressions of these shifts in ( 14.71) we find, in the center-of-mass 
frame, 

44 G 3 m 4 v 2 / r \ 

-^-V-ln - )x<'J>+0(8), (6.5) 



Vu = 

J \' / 

where r' is given by m\nr' Q = m\ In r\ + m 2 In r' T This modification of the quadrupole moment is 
seen to exactly cancel the \n(r/r' Q ) dependence of the mass quadrupole moment in SH coordinates 
as given by (13.1 1 ) — (13 . 21 ) . Thus in the MH gauge the r' Q dependence of the mass quadrupole moment 
vanishes as expected. The rest of the expression of the moment remains exactly the same as in SH 
coordinates, Eqs. 



(13 . 1 1) — (13 .21) . and will not be repeated here. 
Next we must take into account the fact that when computing the third time derivative of the 
quadrupole moment, needed in the expression of the flux, the acceleration in MH coordinates is 
modified. We get 



d 



3vf 



dS^a^ 
dt 



+ 1^2 + (9(8) . 



(6.6) 



where the dots mean the time derivative. The first term is the third time derivative of the direct 
modification of the quadrupole moment, Eqs. (16.4I) - (16.5I) . and the extra terms come from the modi- 
fication of the accelerations which are given by (14.81) . On the other hand, all the other contributions 
coming from the higher multipole moments and their derivatives remain unchanged. We then find 



2C 

s,r = -— 5 [•/,;«>,(•/,;) + a8)|. 



5c 5 



(6.7) 



With the explicit expression of the shifts one finally obtains the modification of the 3PN energy 
flux in the MH coordinates as (thus, Tun = ^sh + S^F) 



5 e T = 



1408 G b mV 



15 



C n r 7 



(6.8) 



[Of course the result agrees with the one we would obtain from using directly Eqs. (16.2l) - (16.3l) .l 
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B. The ADM coordinates 



Many related numerical relativity studies are in ADM (or ADM-type) coordinates and hence 
for future applications we wish to provide the explicit expression for the 3PN energy flux in ADM 
coordinates. To transform the energy flux we require the shift or contact transformation of the 
trajectories connecting the SH coordinates (with log terms) and the ADM coordinates. We recall 
that the ADM and SH coordinate systems agree at 1PN order inclusively, so that the contact 
transformation is composed of 2PN and 3PN terms. Hence the calculation is more involved than 
for the MH coordinates (for which only the modification of the quadrupole moment I t j played a 

role), and we must come back to the general formulas (I6.2| )- (|6.3I) . Note that the remainder O (^ 2 ) 
in Eqs. (I6.2| )- (|6.3I ) is of order 4PN which is still negligible in the transformation to ADM type 
coordinates. 



The relative shift £L linking SH and ADM coordinates, x' 



x' + £' 12 , is given in 114311 as 
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(6.9) 



from which we deduce, applying Eqs. (16.31) . the transformation of variables necessary to compute 
the ADM energy flux: 



5 f r 



Gm 



vH r K Rr ) +q rU +3v 



Gmv f 4 /l 11 \ 2 J 39 99 
+ — T~ \v \ - - —v + r v \ — - + — v 
c 6 1 U 8 / I 16 16 



23 73 
"148" 48 V 



• 4 



Gm 9 /451 3 



+■ 



Gm 



5 \ G 2 m 2 
-v + — — 



Gm 
~r 



r<v 



283 
16 



■t v ) + ' 2 (t v 



2773 22, (r\ 21 , 
+ — In — - — jt 



Gmv 4 / 39 99 

+ —^r\v\ + — v \ + r V 

c 6 r 1 I 8 8 



Gm 
r 

163 
2? 



280 
1 1 

"4 + 2 y 



443 



32 



23 73 

-v\ + r\ + — v 

24 ) I 12 12 



(6.10a) 



For simplicity we use the same notation £' 12 as for the shift between the SH and MH coordinates. 
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(6.10b) 



(6.10c) 



The above equations provide the 3PN generalization of Eq. (4.6) of lfl9ll . They also incorporate 
the corrected transformation between ADM and harmonic coordinates at 2PN, as given in [12511 . 

Using the latter expressions, one finds that the SH energy flux is changed by corrections at 2PN 
and 3PN relative orders given by (using a notation similar to that introduced above - i.e. in which 
the variables r, r 2 and v 2 are considered as dummy variables 7 ) 
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The examination of the coefficient of ln(r/r' Q ), given by the two last terms of (16.1 II) . reveals that this 
coefficient is the same as in the contact transformation from SH to MH, given by (16.81) . Therefore 
the contact transformation from SH to ADM exactly cancels out the logarithms of SH coordi- 
nates, and the final flux in ADM coordinates is free of In r' Q . This is consistent with the general 
understanding that the In r' is a feature of a particular harmonic coordinate system and that ADM 
coordinates do not yield complications associated with such logarithms (the cancellation of the 
In r' {) terms usually provides a useful internal check of the computations). 



VII. THE GENERALIZED QUASI-KEPLERIAN REPRESENTATION 

Before we discuss the calculation of the orbital average of the energy flux in Sect. IVIIIL we 
must summarize the 3PN generalized quasi-Keplerian representation of the binary motion recently 
obtained by Memmesheimer, Gopakumar and Schafer Il37n . Indeed, the main application of the 
present computation is the evolution of the orbital elements under GW radiation reaction to 3PN 
order. This requires one to average over an orbit the instantaneous expressions for the energy flux 
obtained in Sect. |V] Averaging over an orbit is most conveniently accomplished by the use of an 
explicit solution of the equations of motion. The generalized quasi-Keplerian (QK) representation 



of the motion at 3PN order 113711 constitutes an essential input for the computations to follow. 

The QK representation was introduced by Damour and Deruelle 11 1 8M to discuss the problem 
of binary pulsar timing at 1PN order, where relativistic periastron precession first appears and 
complicates the simpler Keplerian picture. This elegant formulation also played an important role 
in our computation of the hereditary terms in Paper I where we provided a summary of it. The 



2PN extension of this work in the ADM coordinates was next given in Refs. Il21l. l22l l23h and we 
shall now use the 3PN parametrization in ADM and MH coordinates [371]. 
The radial motion is given in parametric form as 8 

r =a r (1 - e r cos u) , (7.1a) 
I =u - e t sinw + f t sin V + g t (V - u) + i t sin2V + h t sin3V , (7.1b) 

while the corresponding angular motion is 

4>- (f) P 



K 



V + sin2V + sin3V + sin4V + h# sin5V . (7.2) 



The four angles V, u, £ and (f> are respectively the true anomaly, the eccentric anomaly, the mean 
anomaly and the orbital phase (V, u and I are measured from the periastron, and we denote by 
P the value of (p at periastron). The mean anomaly is proportional to the time elapsed since the 



For convenience in this paper we adapt somewhat the notation with respect to the one in Ref. IB7ll . 
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instant tp of passage at periastron, 

€ = n (f - t P ) , (7.3) 
where n = 2n/P is the mean motion and P is the orbital period. The true anomaly V is given by 



V = 2 arctan 



tan - 



\-eJ 2 



(7.4) 



In the above a r represents the semi-major axis of the orbit, and e r , e t , e<p are three kinds of eccen- 
tricities, labelled after the coordinates t, r and 0, and which differ from each other starting at the 
1PN order. The constant K is linked with the advance of periastron per orbital revolution, and is 
given by K = <D/(27r) where <D is the angle of return to the periastron. The notation k = K - 1 
for the relativistic precession is used in Paper I and will also be useful here. The orbital elements 
ft, ft/,, Qt, 9t/>, ■ ■ ■ parametrize the 2PN and 3PN relativistic corrections, as will be clear from their 
expressions below. [More precisely, f t , fy, g t , g$ are composed of 2PN and 3PN terms, but i t , i^, 
h t , start only at 3PN order.] 

Crucial to the formalism are the explicit formulas for all the orbital elements and all the coeffi- 
cients in Eqs. (17.11 ) above in terms of the 3PN conserved orbital energy E and angular momentum 
J (divided by the binary's reduced mass). Recall that the construction of a generalized quasi- 
Keplerian representation exploits the fact that the radial equation - which is given by Eq. (2.1a) 
in paper I - is a polynomial in 1/r (of seventh degree at 3PN order). Therefore the presence 
of logarithmic terms in the SH coordinates at 3PN order obstructs the construction of the QK 



parametrization (at least by this method) and Ref. 113711 obtained it in coordinates avoiding such 
logarithms, namely the MH and ADM coordinates. In both ADM and MH coordinates the QK 
representation takes the same form given by Eqs. (17.11) - (17.21 ), but of course the equations linking 
the orbital elements to E and J are different. These have been obtained as post-Newtonian series 



up to 3PN order in Ref. [37]. Since they form the basis for our computation of the average energy 
flux, we provide the complete relations here. 

For convenience in the present paper we introduce a PN parameter which is directly linked to 
the energy E and defined by 

2E 

s = -—. (7.5) 

c z 

[Recall that E < for gravitationally bound orbits.] The equations to follow will then appear 
as PN expansions in terms of s = 0(2). Also, we find useful to define, in place of the angular 
momentum J, 

2EJ 2 
(Gmf 

We have j = -2Eh 2 in terms of the more usual definition h = J/(Gm). This parameter is at 
Newtonian order, j = c 2 sh 2 = O(0). The point is now to give all the orbital elements as PN series 
in powers of s with coefficients depending on j (and the dimensionless reduced mass ratio v). In 
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ADM coordinates these are given by Il37ll 
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The latter expressions are specific to the ADM coordinates and we want to give now the corre- 
sponding expressions in MH coordinates. However we recall first an important point related to the 



use of gauge invariant variables in the elliptical orbit case as stressed by Ref. [37]. Indeed Damour 
and Schafer [21] showed that the functional form of n and K = 0>l(2n) as functions of gauge 
invariant variables like s and j are identical in different coordinate systems. Hence the expressions 
in MH coordinates of these two parameters are the same as in ADM coordinates, 

n ee n MB =n ADM , (7.8a) 
K = K MU =K ADM . (7.8b) 



This prompted Ref. H37I1 to suggest the use of n and k = K - 1 as two gauge invariant variables in 
the general orbit case. 9 In the present work we propose to use a variant of the former variables. 
Namely, instead of working with the mean motion n we shall systematically use the orbital fre- 
quency a> = Kn as defined in a general context in Sect. IIA of Paper I, and define as a gauge 
invariant post-Newtonian parameter 

. C7.9, 

This choice constitutes the obvious generalization of the gauge invariant variable x used in the 
circular orbit case and will thus facilitate the straightforward reading out and check of the circular 
orbit limit. The parameter x is related to the energy and angular momentum variables s and j up 
to 3PN order by 

f 1 / 5 1 ol 

x = e{\ +e -- + — v + 2- 
l \ 4 12 j 
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(7.10) 



The other orbital elements are not gauge invariant and therefore their expressions in MH coordi- 
nates differ at 2PN and 3PN orders from those in ADM coordinates. We conclude by giving here 



Actually Ref. |37|] used x M gs = (Gmw/c 3 ) 2/3 together with k' = k/3. 
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all the needed differences B7I1 . 
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iT - C DM =-^(-1 + 7) 2 r 3 v(-149 + 198v) . (7.111) 

Finally we note that in the case of a circular orbit the angular momentum variable, say j @ , is related 
to the constant of energy s by the 3PN gauge-invariant expansion 

j Q = 1 + - (9 + v) + — (81 - 32v + v 2 ) + 7^(2835 - 7699v + 246vtt 2 + 96v 2 + 3v 3 ) , (7.12) 
4 16\ /192V / 

which is easily deduced using either MH or ADM coordinates. This expression can be used to 
compute all the orbital elements for circular orbits and we can check that all of the eccentricities 
e r , e t or e* are zero. 



VIII. ORBITAL AVERAGE OF THE 3PN ENERGY FLUX 

To average the energy flux over an orbit we will require the use of the previous 3PN quasi- 
Keplerian representation of the motion. Consequently, the averaging is only possible in MH or 
ADM coordinates without the logarithms as discussed before. The average of the (instantaneous 
part of the) energy flux is defined by 

1 C P 1 C 2n dl 

<^ nst > = - dtTi nst = — \ du — T mA . (8.1) 
P Jo 2 7r Jo du 

As we have seen the energy flux (12.61) is made of instantaneous terms and hereditary (tail) terms. 
The hereditary terms have already been computed and averaged in Paper I. 

Using the QK representation of the orbit discussed in Sect. IVIH we can re-express the energy 
flux ^ nst [or, more exactly, {dt/du) TinstL which is a function of its natural variables r, r and v 2 , as a 
function of the frequency-related parameter x defined by Eqs. (17.9D - ti7.10D , the "time" eccentricity 
e t and the eccentric anomaly u. w We note that in the expression of the energy flux at the 3PN 
order there are some logarithmic terms of the type ln(r/r ) even in MH coordinates. Indeed, we 



10 Ref. IU9I1 uses Gm/a r and e r while 02511 employs Gmn/c 3 and e,. We propose the use of x = (Gmco/c 3 ) 2 ^ for 
reasons outlined in the previous Section. The choice of e, rather than say e r is a matter of convenience since it 
appears in the Kepler equation which is directly dealt with when averaging over an orbit. 
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recall that the MH coordinates permit the removal of the log-terms ln(r/r ), where r' Q is the scale 
associated with Hadamard's self-field regularization, but there are still the terms ln(r/r ) which 
involve the constant r entering the definition of the multipole moments for general sources. As a 
result we find that the general structure of Tm A (in MH or ADM coordinates) is 
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(8.2) 



where the coefficients a^, 0n, Jn so defined are straightforwardly computed using the QK 
parametrization (they are too long to be listed here). It is worth noting that the yS^'s correspond to 
all the 2.5PN terms while the y w 's represent the logarithmic terms at order 3PN. The dependence 
on the constant In r has been included into the coefficients a N 's. To compute the average we have 
at our disposal some integration formulas. First of all, 
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which shows that in the final result there will be no terms (of the "instantaneous" type) at 2.5PN 
order. The 2.5PN instantaneous contribution is proportional to r and vanishes after averaging since 
it includes only odd functions of u. Next, we have 
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which can also be formulated with the help of the standard Legendre polynomial P^-x as 
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Finally for the log-terms we have a less trivial formula but which takes a structure similar as in 
Eq. (I8ll) . namely 
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(8.7) 



A. Orbital average in MH coordinates 

The expression for the instantaneous energy flux in MH coordinates is given by Eqs. (15.11) - 
(15.21 ) together with the modification (16.81 ) for transforming to MH coordinates. Implementing all 
the above integrations, the flux can be averaged over an orbit to order 3PN extending the results 
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of Iil9l1 at 2PN 11 . The result is presented in the form 
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where the "instantaneous" post-Newtonian pieces depend on v and the time eccentricity e t in 
MH coordinates (note that e t = ef H here), and read 12 
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(8.9b) 



(8.9c) 



11 Results of IU9I1 are given in ADM coordinates. 

12 The Newtonian coefficient is nothing but the Peters & Mathews [ 1] "enhancement" function of eccentricity 
f(e t ) = (1 + ||e 2 + ||e^)/(l - e 2 ) 7/2 , called that way because it enhances the numerical value of the orbital decay 
of the binary pulsar by gravitational radiation (viz the orbital P). 
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For ease of presentation we have not put a label on e t to indicate that it is the time eccentricity 
in MH coordinates. Of course, since x is gauge invariant, no such label is required on it. It 
is important to keep track of this fact when comparing formulas in different gauges, as we will 
eventually do. 

The last term in the 3PN coefficient given by Eq. ( |8.9d|) is proportional to some logarithm 
which directly arises from the integration formula J8.6lHl8.7l) . Inside the logarithm we posed 
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exhibiting the dependence of the instantaneous part of the 3PN energy flux upon the arbitrary 
constant length scale r . Only after computing the complete energy flux can one discuss the 
structure of the logarithmic term in the energy flux and the required cancellation of the lnr . 
Therefore we now add the hereditary contribution to the 3PN flux which has been computed in 
Paper I. From Eq. (6.2) in Paper I we write the result as 
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where the "hereditary" post-Newtonian coefficients (starting at 1 .5PN order) read 
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The function F(e t ) in factor of the logarithm in the 3PN coefficient does admit a closed analytic 
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form which was determined in paper I as 
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On the other hand paper I found that the four "enhancement" functions of eccentricity (p(e t ), iff(e t ), 
£(e t ) and K(e t ) very likely do not admit any analytic closed-form expressions. Numerical plots of 
the four enhancement factors <p(e t ), \]/(e t ), 6{e t ) and K(e t ) as functions of eccentricity e t have been 
presented in Paper I. The coefficients in Eqs. (18.121) have been introduced in such a way that the 
circular orbit limit of all the functions F(e t ) and <p(e t ), • • • , K{e t ) is one. 

Finally, the PN coefficients in the total averaged energy flux T in MH coordinates are given by 
the sum of the instantaneous and hereditary contributions, say 



tMH 



''nPN 



(8.14) 



We notice that up to 2.5PN order there is a clean separation between the instantaneous terms 
which are at even PN orders (recall that there is no 2.5PN term in the averaged flux), and the 
hereditary terms which appear at odd PN orders and are specifically due to tails (i.e. , H 1 M 5 p N and 
< H^ 1 5 p N ). On the contrary, at 3PN order - and, indeed, at any higher PN order - there is a mixture 
of instantaneous and hereditary terms. The 3PN hereditary term *H^S is due to the so-called GW 
tails of tails (see Paper I). 

The analytical result (18.131) is crucial for checking that the arbitrary constant xq disappears 
from the final result, namely from the 3PN coefficient 9C|j£. Indeed we immediately verify from 
comparing the last term in Eq. (|8.9d| ) with Eq. (I8.12cl) and the explicit expression (18.131) of F(e t ) 
that x cancels out from the sum of the instantaneous and hereditary contributions, extending to 
non circular orbits this fact which was already observed for the circular case in Ref. lf30h . Finally 
the complete 3PN coefficient (independent of xq) reads 
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The 1.5PN and 2.5PN coefficients are only due to tails, thus 

<K^=Anip(e t ), 
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The Newtonian, 1PN and 2PN coefficients reduce to their instantaneous contributions I^ H , 
and already given in Eqs. (18.91) . 

Since the enhancement functions tp(e t ), ifr(e t ), £(e t ) and K{e t ) reduce to one in the circular case, 
when e t = 0, the circular-orbit limit of the energy flux is immediately deduced from inspection of 
Eqs. dlill) and (f8T6l) as 
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This limiting case is in exact agreement with Eq. (12.9) of 113011 (after taking into account the values 
of the ambiguity paramete 
the flux in the circular-orb 
becomes gauge invariant. 



of the ambiguity parameters A = -3^5 and 6 = --^jf computed in Refs. ||32Ll33L l34ll N ). Notice that 
the flux in the circular-orbit limit (18.171) depends only on the parameter x and hence its expression 



B. Orbital average in ADM coordinates 

We start from the expression for the instantaneous energy flux in ADM coordinates as given 
by (16.1 II) . employ the appropriate 3PN QK representation and follow the procedure for performing 
the average as outlined in the previous Section. We find that the fi N 's, and y^'s in ADM coordinates 
[cf. Eq. (18.21) 1 are exactly the same as in MH coordinates; the c/v's, however, are different in general 
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(except for an). The result for the average energy flux in ADM coordinates is of the form 
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where the coefficients depend on the time eccentricity in ADM coordinates (hence e t = ey here) 
and on v, and read 

1 



J 



ADM 
N 



I 



ADM 
1PN 



1 



73 



37 



1 + — e; + — e, 
24 ' 96 ' 

1247 35 2 / 

-^36--l2 V + 



10475 1081 



I 



10043 311 



f ADM 
x 2PN 



\ 384 
1 



y\ + e: 



12 
203471 
9072 



f \ 672 
2179 851 



36 



+ 



1792 576 
12799 65 



v H y 

504 18 



3866543 4691 5935 



+ 



18144 2016 



■v + 



369751 3039083 



24192 8064 
1302443 215077 



■y + 



54 
247805 
864 
185305 . 



' \ 16128 1344 ' 1728 
86567 9769 21275 



+ 



4 

'\64512 4608 6912 

r :[35 ^ 2 /6425 1285 \ 

Jl-e 2 — -7y + < — — y 

> ' [ 2 ' \ 48 24 / 



5065 1013 



v +e 



J 



ADM 
3PN 



64 32 

[ 2193295679 
I 9979200 



185 37 



(1 - e 2)13/2 



+ 



96 48 , 
8009293 41tt 2 



54432 



64 



y - 



209063 
3024 



y 2 - 



775 
324 1 



+et 



2912411147 



2851200 



+ 



249108317 31255 



-7T 



+e;\- 



4520777971 



+e 



+ 



13305600 
/ 3630046753 
\ 26611200 
/ 21293656301 
' \ 141926400 



+ 



108864 1536 
473750339 7459^ 2 



y - 



3525469 , 53696 
v - 



108864 1024 
8775247 78285;r 2 



y + 



6048 
697997 



+ 



145152 4096 
36646949 4059tt 2 



y + 



129024 



4096 



v + 



576 
31147213 
12096 
85830865 



y 2 - 



243 
10816087 
7776 



y 2 - 



983251 
648 



193536 



v 2 - 



4586539 
15552 1 



(8.19a) 



(8.19b) 



(8.19c) 



8977637 9287 



8977 



+ 



v + 



11354112 48384 55296 



y 2 - 



567617 



124416 



+ 



14047483 
151200 



+ 



165761 287tt 2 



1008 



192 



y + 



455 



32 



+ e 



7 36863231 
' 100800 



+ 



14935421 52685?r 2 



+ 4 



759524951 



+ 



+ e 



403200 
/ 139966120 3 
' \ 2419200 
185 1073 



+ 



6048 
31082483 

8064 4 
40922933 



48384 



+ e: 



+ 



48 
1712 



105 



+ 



288 

1 4552 
~63~ 



-v + 



407 



e 2 t + 



288 

553297 4 
1260 Sf 



4608 
41533tt 2 
6144 

15177T 2 

f 9216 



v + 



43559 

IT 1 



V + 



V + 



3 03985 

288 

73357 , 

Y 

288 



187357 
1260 



10593 c 

c? 

2240 ' 



In 



x 1 + ^J\- e 
~o 2(1-^) 



(8.19d) 



We recall that the Newtonian and 1PN orders are the same in MH and ADM coordinates (the 
coefficients i"^ DM and i"f™ agree with their MH counterparts). On the other hand, adding up the 
hereditary contribution (18.1 1|) - (I8. 121 ) [which is the same in MH and ADM coordinates], we obtain 
the total 3PN coefficient < K^ A , ana l°g° us to Eq. (18.151 ) but in ADM coordinates, 
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in which again e t = e^ DM . A useful internal consistency check of the algebraic correctness of dif- 
ferent coordinate representations of the energy flux, is the verification that the equality of Eqs. (18.91) 
and (18.191 ) holds if and only if we have the transformation between the time eccentricities and 
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(There is no ambiguity in not having a label on the e t in the 2PN and 3PN terms above.) We 
find that the relation (18.211) is perfectly equivalent to what is predicted from using different QK 
representations of the motion, namely Eq. (17.1 lcl) together with (17.101) . 



C. Gauge invariant formulation 

In the previous section, the averaged energy flux was represented using x - a gauge invariant 
variable defined by (17.91 ) - and the eccentricity e t which however is coordinate dependent (but is 
useful in extracting the circular limit of the result). In the present Section we provide a gauge 
invariant formulation of the energy flux. 

Perhaps the most natural choice is to express the result in terms of the conserved energy E 
and angular momentum J (per unit of reduced mass), or, rather, in terms of the pair of rescaled 
variables (s, j) defined by Eqs. (17.51) and (I7.6I ). However there are other possible choices for 
a couple of gauge invariant quantities. As we have seen in Eqs. (17.81) the mean motion n and 
the periastron precession K are gauge invariant so we may define as our first choice the pair of 
variables (x, i), where we recall that x is related to the orbital frequency a) = Kn by Eq. (17.91 ), and 
where we define 

(8.22) 

k 

with k = K - 1. Here we have introduced a factor 3 so that i reduces to j in first approximation 
{i.e. when s —> 0). To 3PN order this parameter is related to the energy and angular momentum 
variables e and j by 
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We have performed two calculations of the gauge-invariant result, in terms of the variables (x, 
i), starting from the expression of the averaged flux in either MH and ADM coordinates. The 
instantaneous part of the flux takes the form 
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in which the PN coefficients are polynomials of t and the mass ratio v and given by 
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Similarly, we can also obtain the equivalent expression of the flux in terms of the rescaled variables 
(e, j) defined by Eqs. £73]) and (171)1) . 

The hereditary part of the flux given by (18.1 1I) - (I8. 121) is straightforwardly added. In this part 
we have simply to replace e t by its expression in terms of x and l at the 1PN order, namely 
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(At this order there is no difference between MH and ADM coordinates.) Note also that with the 
latter choice of gauge-invariant variables the circular-orbit limit is not directly readable from the 
expressions. However, it can be easily obtained by using the expression for the variable j Q as 
reduced to circular orbits in terms of e, Eq. (17.121) . 



IX. THE TEST PARTICLE LIMIT OF THE 3PN ENERGY FLUX 

An important check on our result is the test particle limit for which the energy flux in the 
eccentric orbit case is available (to second order in the eccentricity) from computations based 
on perturbation theory around a Schwarzschild background. We compare the end result of our 
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computation - composed of the instantaneous terms and the hereditary terms computed in paper I 
- with the result obtained in Ref. lf4vh . Thus, we take the test particle limit of our result (i.e. 
v = \x I m — » 0), say in the form given by Eqs. (18.8lHl8.9l ) in which e, = ef 111 , and expand it in 
powers of e t retaining only terms up to e 2 . The instantaneous contribution to the energy flux in the 
test mass limit is then given by 
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On the other hand, the hereditary contribution has been reported in Eqs. (18.1 1|) -( 18. 121 ) and admits 
the test-mass limit 
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To proceed further, all the enhancement functions should be expanded up to power e 2 . This is easy 
for F(e t ) which is known analytically from Eq. (18.131) and we have 
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F(e t ) = l + ^e 2 + 0(4). 



(9.3) 



The other enhancement functions are only known numerically for general eccentricity. We have 
however succeeded in obtaining analytically their leading correction term e 2 by implementing 
our calculation of the tails in paper I at order e 2 from the start. The results we thereby obtained 
[Eqs. (6.8) of paper I] are 
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[We do not need g(e t ) here since it is in factor of a v-dependent term.] Our final result to 0(v) and 
0(e^ is therefore 
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x 3 j + <9(ef) + (9(y)|. (9.5) 



The above expression is in terms of our chosen eccentricity e t . One should note that the 
"Schwarzschild" eccentricity e appearing in the black-hole perturbation theory Il47ll is a priori 
different from e t ; therefore the above result can only be compared modulo a transformation of 
these eccentricities. We find that indeed Eq. (19.51 ) is equivalent to the black-hole perturbation 
result given by Eq. (180) of |47|], if and only if the two eccentricities are linked together by 



?f = e 2 (l - 6x + Ax 2 - 8jc 3 ) 



(9.6) 



(Recall that e, = e, here.) 



X. CONCLUDING REMARKS 

The instantaneous contributions to the 3PN gravitational wave luminosity from the inspiral 
phase of a binary system of compact objects moving in an elliptical orbit is computed using the 
Multipolar post-Minkowskian wave generation formalism 13 . The non-trivial inputs for this calcu- 
lation include the mass octupole and current quadrupole at 2PN order for general orbits and the 
3PN accurate mass quadrupole. Using the 3PN quasi-Keplerian representation of elliptical orbits 
obtained recently the flux is averaged over the binary's orbit. The instantaneous part of the energy 
flux is computed in the standard harmonic coordinate system (with logarithms). For technical rea- 
sons the average over an orbit of the instantaneous contributions is presented in other coordinate 
systems: Modified harmonic coordinates (without logarithms) and ADM coordinates. Alternative 
gauge invariant expressions are also provided. Supplementing the instantaneous contributions 
of this paper by the important hereditary contributions arising from tails, tails-of-tails and tails 
squared terms calculated in paper 1 1|3^1 . the complete energy flux has been obtained. 



For binaries moving on circular orbits the 3PN energy flux agrees with that computed in [13011 . 
However the circular-orbit results are known to the higher 3.5PN order lL30h . The extension of the 
3.5PN term to eccentric orbits would be interesting, but some uncomputed modules remain in the 
general formalism to compute the multipole moments for general sources required for the 3.5PN 
generation in the eccentric orbit case. We leave this to a future investigation. 
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The instantaneous part of the 3PN gravitational wave flux of angular momentum and linear momentum from inspi 



railing compact binaries moving on elliptical orbits has been computed 11481 14911 . 
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to the results of this paper are done with the software MATHEMATICA. 
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